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The transient behavior of an ecosystem with N random interacting species in the presence 
of a multiplicative noise is analyzed. The multiplicative noise mimics the interaction with 
the environment. We investigate different asymptotic dynamical regimes and the role of the 
external noise on the probability distribution of the local field. 

§1. Introduction 

Population dynamics attracted a lot of attention in recent years and became 
the object of many studies as well by biologists as by physicists Tools 
developed in the context of nonequilibrium statistical physics to analyze nonequi- 
librium nonlinear physical systems provide new insights and at the same time new 
approaches to study biological systems. Biological population dynamics has many 
interesting, and still not solved, problems such as pattern formationP^^^ the 
role of the noise on complex ecosystem behaviour, and the noise-induced effects, 
such as stochastic resonance, noise delayed extinction, quasi periodic oscillations 
etc.lBajmOOOJXSIJXSIJniJXSI The dynamical behavior of ecological systems of 
interacting species evolves towards the equilibrium states through the long, slow, and 
complex process of nonlinear relaxation, which is strongly dependent on the random 
interaction between the species, the initial conditions and the random interaction 
with environment. A good mathematical model to analyze the dynamics of iV bio- 
logical species with spatially homogeneous densities is the generalized Lotka-Volterra 
system with a Malthus-Verhulst modelization of the self regulation mechanism, in 
the presence of a multiplicative noiseP*'*=*'*=^' By neglecting the fluctuations of the 
local field we derive a quasi-stationary probability of the populations. We obtain 
the asymptotic analytical expressions for different nonlinear relaxation regimes, and 
we analyze the role of the multiplicative noise on the probability distribution of the 
local field. 

§2. The model and results 

The dynamical evolution of our ecosystem composed by N interacting species 
in a noisy environment and in the presence of an absorbing barrier is described by 
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the following Ito stochastic differential equation 



dxi(t) 



9i{xi{t)) + ^2 J ij x j(t) ) dt + V^dwi 



Xi (t), i = l,...,N (1) 



where Xi(t) > is the population density of the i th species at time t and the function 

9i{xi{t)) 

g i (x i (t))=(a + ^-x i (t) (2) 

describes the development of the i th species without interacting with other species. 
In Equation a is the growth parameter, the interaction matrix Jij models 
the interaction between different species (i ^ j), and Wi is the Wiener process 
whose increment dwi satisfies the usual statistical properties < dwi(t) > = 0, and 
< dwi(t)dwj(t') >= 5ij5(t — t')dt. We consider a random asymmetric interaction 
matrix J^, whose elements are independently distributed according to a Gaussian 
distribution with (J^) = 0, (JijJji) = 0, and a 2 = J 2 /N. The term JijXiXj is the 
loss or the growth rate of species % due to interaction with species j, when Jjj < or 
Jij > respectively. With this choice of interaction matrix our ecosystem contains 
50% of prey-predator interactions (Jij < and Jij > 0), 25% competitive interac- 
tions (Jij < and Jij < 0), and 25% symbiotic interactions (Jij < and Jij > 0). 
We consider all species equivalent so that the characteristic parameters of the ecosys- 
tem are independent of the species. The random interaction with the environment 
(climate, disease, etc..) is taken into account by introducing a multiplicative noise 
in the Equation (0). The solution of the dynamical Equation Q is given by 

/ \ Xi(0)zi(t) , . 

xAt) = n ' P ' , 3) 

where 

z At) = e^+V^WiW+Jo d,tfrn(t') 

and 



Vi(t) = ^JijXj(t) (5) 

is the local field acting on the i th population and represents the influence of other 
species on the differential growth rate. We note that the dynamical behavior of the 
i th population depends on the time integral of the process Zi(t) and the time integral 
of the local field. For a large number of interacting species we can assume that 
the local field rn(t) is Gaussian with zero mean and variance a 2 . = (rjj) = J 2 (xf). 
As a consequence, in the absence of external multiplicative noise, from the fixed- 
point equation Xi(a — Xi + rji) = 0, the stationary probability distribution of the 
populations is the sum of a truncated Gaussian distribution at Xj = (xi > 



Transient population dynamics 



3 



always) and a delta function for extinct species. The initial values of the populations 
Xi(0) have also Gaussian distribution with mean value (a^(0)) = 1, and variance 
0-2 = 0.01,0.03,0.05. 

The interaction strength between the species J determines two different dy- 
namical behaviors of the ecosystem. Above a critical value J c = 1.1, the system 
is unstable and at least one of the populations diverges. Below J c the system is 
stable and asymptotically reaches an equilibrium state. The equilibrium values of 
the populations depend both on their initial values and on the interaction matrix. 
If we consider a quenched random interaction matrix, the ecosystem has a great 
number of equilibrium configurations, each one with its attraction basin. For an 
interaction strength J = 1 and an intrinsic growth parameter a = 1 we obtain: 
< Xi >= 1.4387, < x\ >= 4.514, and a 2 x . = 2.44. These values agree with that ob- 
tained from numerical simulation of Equation (JJ) . From the Fokker-Planck equation 
associated to the Langevin equation (JJ) 



d_ 

dt 



P(Xi,t) 



dxi 



e_0_ 
2dx~ 



(a + - - Xi + m)xi 



P( Xi ,t) 



(6) 



we obtain a quasi stationary distribution by neglecting the fluctuations of the local 
field in the asymptotic regime 



dP{ Xl ) 

P(Xi) 



2 r 

e 



a + rji 



Xi 2 



dxi 



(7) 



that is 



P(xi) = N Xi( p)exp 
where the normalization factor is 



- ((a +r)i - |) ba.Xi - x, 



N 



Xi(0) 



P Xi(,0) e 



2^(0) 



O(Xi) 



(8) 



0) 



x .( ) 7 («+^-f) 

and is the Heaviside unit step function. Now we focus on the statistical properties 
of the time integral of the i th population Xi{t) 



X l (t)= / dt' Xi (t'), 



(10) 



in the asymptotic regime. From Eqs. © and (|10|) we have 



Xi(t) = In 



1 + Xi{0) [ dt'exp[at' + ^Wi(t') +rn(t')] 
Jo 



ID 



We use the same approximation of the mean field interactionf^' and after differen- 
tiating Eq. (JTT|) , we get the asymptotic solution of X{ (t) as 



Xi(t) ~ln 



Xi (o) e ^ma a:j (t)+w j (*) / dt > e at' 

Jo 



(12) 
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where w maXi (t) = sup 0< t> <t w(t') and rj maXi (t) = sup 0<t i <t ri(t')- Eq. (H2J) is valid for 
a > 0, that is, when the system relaxes toward an equilibrium population and at the 
critical point. After making the ensemble average, we obtain for the time average of 
the i th population Xi 

, a = 0, (13) 



and 

where N w and N„ are variables with a semi-Gaussian distribution^ and N„ must 
be determined self-consistently from the variance of the local field (Eq. ©). We 
obtain, consistently with mean field approximation, the typical long time tail be- 
havior (t -1 / 2 ) dependence, which characterizes nonlinear relaxation regimes when 
a > 0. When the system relaxes toward the absorbing barrier (a < 0), the time 
average of the i th population (-Xj) is a functional of the local field and the Wiener 
process. We have also analyzed the dynamics of the ecosystem when one species is 
absent. Specifically, we considered the cavity field, which is the field acting on the 
i th population when this population is absent. The probability distributions for both 
local and cavity field have been obtained by simulations for a time t = 200 (a. u.) 
in absence of external noise, and for different species. We found that the probability 
distributions of the cavity fields differ substantially from those of local fields for the 
same species, while in the presence of noise the two fields overlap (see Figs. 1 and 
2). 



N w Vet + lni + (In Mo)]) 



N w Vei + (a + N v )t + (in 



Xi(0) 



a 



, a > 0, 



(14) 




This overlap is different for different species and depends on the variance of 
the initial species distribution. This strange behavior, found for some populations, 
is reminiscent of the phase-transition phenomenon, and it is related to the follow- 
ing peculiarities of our dynamical system: (i) all the populations are positive; (ii) 
different initial conditions drive the ecosystem into different attraction basins; and 
(iii)the complex structure of the attraction basins of our dynamical system. While, 
in the presence of noise, all the populations seem to be equivalent from the dynamical 
point of view, some populations, in the absence of external noise, have an asymptot- 
ical dynamical behavior such that they significantly influence the dynamics of other 



Transient population dynamics 



5 



species. 



§3. Conclusions 



We analyzed the nonlinear relaxation of an ecosystem composed by N interacting 
species. We obtain the quasi-stationary probability distribution of the population 
in the presence of multiplicative noise. By using an approximation of the integral 
equation, which gives the stochastic evolution of the system, we obtained asymp- 
totic behavior for different nonlinear relaxation regimes. We observe an interesting 
phenomenon: the local and the cavity fields, whose probability distributions are 
different in the absence of noise, coincide for some populations in the presence of 
multiplicative noise. This phenomenon can be ascribed to the peculiarity of the 
dynamical system (xj > 0, always), the influence of different initial conditions on 
the asymptotic regime, and the complex structure of the attraction basins of our 
ecosystem. 
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